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Abstract 

We present an exact solution of the deterministic Kardar-Parisi-Zhang (KPZ) 
equation under the influence of a local driving force /. For substrate dimen- 
sion d < 2 we recover the well-known result that for arbitrarily small / > 0, 
the interface develops a non-zero velocity v{f). Novel behaviour is found in the 
strong-coupling regime for d > 2, in which / must exceed a critical force fc in 
order to drive the interface with constant velocity. We find v{f) ^ {f — fc)"'^'^^ 
for / \ fc- In particular, the exponent a{d) = 2/{d — 2) for 2 < < 4, but 
saturates at a{d) = 1 for d > 4, indicating that for this simple problem, there 
exists a finite upper critical dimension du = 4. For d > 2 the surface distor- 
tion caused by the applied force scales logarithmically with distance within 
a critical radius Rc ^ {f — fc)~'^^'^\ where ^(d) = a{d)/2. Connections be- 
tween these results, and the critical properties of the weak/strong-coupling 
transition in the noisy KPZ equation are pursued. 
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I. INTRODUCTION 



This article is concerned with a problem which arises in the study of the Kardar-Parisi- 
Zhang (KPZ) equation of interface growth Namely, we consider the deterministic 
version of the KPZ equation under the influence of a local, constant force. Although this 
problem is in principle rather simple, a detailed investigation has not been carried out 
to our knowledge. The results, although interesting in their own right, assume greater 
significance when related to the weak/strong-coupling transition in the noisy KPZ equation. 
Before entering into further details of this problem, we shall give a brief overview of the 
key properties of the KPZ equation. In, by now, standard notation this equation for the 
evolution of the interface height /i(x, t) (above a (i-dimensional substrate) takes the form 

dth = uV^h + {X/2)(yhf + rj (1) 

where ri{'x,t) is a white noise, generally taken to be gaussian distributed. The Cole-Hopf 
transformation /i(x, t) = (2z//A) ln(w(x, t)) yields a linear equation for w albeit with multi- 
plicative noise: 

dtw = vW'^w + (A/2z/)wr/. (2) 

This equation may be re-expressed by means of the Feynman-Kac formula, such that w 
represents the restricted partition function of a directed polymer in the presence of quenched 
random defects. Naturally, the above equation may also be interpreted as an imaginary-time 
Schrodinger equation for a quantum mechanical (QM) particle in a time-dependent random 
potential {—rf). 

It is well-known [0] that d = 2 plays the role of a lower critical dimension. For d <2, the 
interface will be asymptotically (i.e. for large times) rough for arbitrarily small A. For d > 2 
there exists a transition between asymptotically smooth and rough phases, depending on the 
value of A. To access the rough (strong-coupling) phase, A must exceed a critical value Ac. 
Amongst the most important questions in this field are: i) what are the scaling exponents 
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associated with the strong-couphng phase; ii) does there exist an upper critical dimension du, 
above which these exponents take on mean-field values; iii) what are the critical properties 
of the weak/strong-coupling transition? There are no answers to the first question, bar 
predictions from extensive numerical work (except for the case d = 1 where exact analytic 
results are available 0). As regards the second question, the community seems split between 
beliefs in (i„ = 4 and du = oo. The former belief is based primarily on the predictions 
of various self-consistent 'mode-coupling' theories (see for example @]), whilst the latter 
belief arises from numerical evidence for d-dependent strong-coupling exponents |^ (for all 
accessible d~ 7). With regard to the final question, a recent two-loop renormalization group 
calculation [Q (formulated as an e-expansion, with e = d — 2) predicts that at the transition, 
the roughness exponent x = -|- O(e^) (to be interpreted as logarithmic roughness), whilst 
the correlation length exponent has the form z/ = (1/e) -|- O(e^). For the one-loop versions 
of these results, see Ref. 0. 

We now concentrate on a much simpler problem. As mentioned before, we retain the 
essential non-linearity of the KPZ equation, but replace the noise ?7 by a local applied force: 
?7(x, t) = /o5'^(x). In directed polymer language this corresponds to replacing the bulk 
disorder by a single constant energy columnar defect. In terms of the QM analogue, we have 
replaced the random potential by a delta-function potential with amplitude (— /o)- There 
exists a well-known result 0, which is most easily stated in QM language: for d < 2 the 
particle will have a bound state for arbitrarily small /o > 0, whereas for d > 2 there only 
exists a bound state for /o > /o,c > 0. In terms of the interface, these results imply that 
for d <2 an arbitrarily small force is sufficient to induce a non-zero velocity to the surface, 
whereas for d > 2 the driving force must exceed a critical value in order to achieve this. 
The identification of a lower critical dimension (along with the properties of the system for 
d < 2) may be obtained from a perturbative analysis (in terms of A, or equivalently /o). 

In order to explore the properties of the 'binding/unbinding' transition for > 2, a 
non-perturbative analysis is required - such an analysis constitutes the body of the present 
work. The main results which arise are the following (which we shall describe in terms of 
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the interface formulation). Whereas for d < 2 an arbitrarily small driving force will induce 
a velocity f (/o) to the interface, for d > 2 the applied force must exceed a critical value in 
order to achieve this. The dependence of the velocity on the excess force Sfo = /o — /o,c is 
found to be ~ (5/o)°^'^^ where a{d) = 2/(d - 2) for 2 < d < 4, but saturates at a{d) = 1 
for d > 4. For d < 2 the applied force will create a non-hnear distortion of the surface over 
a critical radius Rc ~ fo^^^'^~'^\ but the amplitude of the distortion is of order unity. For 
d > 2 the distortion scales logarithmically with distance and exists within a critical radius 
Rc ~ {Sfo)~''^'^\ where u{d) = a{d)/2. Therefore, for d > 2 the distortion IS a macroscopic 
rupture of the surface with a logarithmic profile (by 'rupture' we mean a gross distortion of 
the interface, leaving the continuous properties of the surface in tact). We shall leave the 
interpretation of these results for the final sections of the paper. 

The outline of the paper is as follows. In the next section we shall give a precise defi- 
nition of the model to be solved, along with some simple steps which lead to an accessible 
strong-coupling (i.e. non-perturbative) analysis of the model. In section 3, we proceed with 
this analysis and derive the velocity-force characteristics. The spatial profiles of the sur- 
face distortion for various dimensions are evaluated in section 4. Connections between the 
properties of this simple model, and the weak/strong-coupling transition in the noisy KPZ 
model are pursued in section 5. We end the paper with our conclusions, and a description 
of possible extensions to the present work. 

II. DEFINITION AND REFORMULATION OF THE MODEL 

Replacing the stochastic source of the noisy KPZ equation with a local applied force 
leads us to consider the equation 

dth^vV^h+{\/2){Vhf + U5\^). (3) 

As for the noisy problem, we may apply the transformation h{x,t) — {2u/X)\n{w{x,t)) 
which yields the linear equation 
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dtw = uV^w + (A/o/2z/) w 5'^(x). (4) 

Since we are primarily interested in the interface application of this problem, we choose the 
initial condition ti?(x, 0) = 1, which corresponds to an initially fiat interface {h = 0). [The 
natural initial condition for the directed polymer application is w(x, 0) = 6'^{x). We shall 
not consider this problem here.] 

We may formally integrate Eq.(^) using the Green function of the diffusion equation, 
(7(x, t) = (47rz/t)"°'/^ exp(— x^/4z/t), leading to 

t 

w(x, t) = 1 + (A/o/2z/) J dt' g{x, t - t')w{0, t'). (5) 



Therefore, knowledge of the function iy(0,t) is sufficient to determine w throughout all 
space. For convenience we shall define (f){t) = w(0,t). Setting x = in the above equation 

leads to a Volterra integral equation for : 

t 

(Pit) = 1 + / J dt' a{t- t')(f){t') (6) 



where a{t) = and we have defined an effective force / = (A/o/2z/)(47rz/)^'^/^. [Such a 
Volterra equation may also be derived for the more complicated situation of a time-dependent 
force fit). The case where f{t) is a stochastic variable is of particular interest in the directed 
polymer picture ||^.] 

Clearly this equation is ill-defined for d > 2, since the kernel is non-integrable over the 
range (0,t). To regularize the equation we have many possibilities. The most convenient is 
to smear the delta-function force over a small region of size l'^ with a gaussian envelope, i.e. 
we replace /o5'^(x) by /o(7r/^)^'^/^ exp(— x^//^). The length scale / constitutes the smallest 
physical scale in the problem. Assuming that w(x, t) varies slowly over the core region 
|x| < / allows us to rederive the above Volterra equation, where now a{t) = (t + p)-^/'^ 
(after rescaling 1^1 = Z/(4z/)). 

Laplace transforming Eq. (^) and using the convolution theorem yields an explicit solution 
for this Laplace transform £[0(i(:)]. Inverting the transform gives the result 
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ds ^st/P 

T [1 - P-dfD{s)] 



(7) 



where 7 is the usual integration contour along the line s = c + iy where c is a real constant 
chosen to lie to the right of any singularities in the complex plane. The function D{s) is 
given by 



(which may be expressed in terms of the incomplete Gamma function 0). 

We are now in a position to obtain all physical properties of the model from the analysis 
of Eq.(l^). We stress once again that this analysis is non-perturbative in /q (or equivalently 
A). All the results for d < 2 may be obtained by perturbation theory (although one must 
evaluate all orders), whereas none of the results concerning the binding/ unbinding transition 
may be obtained by such a technique for d > 2. 

In the next section we shall concentrate on evaluating the asymptotic (large time) form 
of (f){t) which will enable us to determine whether the interface has a non-zero velocity. The 
spatial profile of the interface distortion requires the evaluation of the convolution of the 
diffusion equation Green function with the function 0(t) (as illustrated in Eq. (^) - this 
analysis is described in section 4. 



The singularities in the integrand of Eq. (^ consist of a branch cut along the negative real 
axis, and (depending on the value of /) a pole in the right-half plane (RHP). The large-time 
behaviour of 0(t) will be dominated by the contribution from the pole (if it exists). One may 
easily convince oneself that the existence of this pole is required for the interface at the forcing 
site to have a non-zero velocity, i.e. for H{t) ~ vt (where H{t) = {X/2u)h{0,t) = In (f){t).) 

By studying the small-s behaviour of D{s) one may show the following properties of the 
integrand in Eq.(0). For d < 2 there always exists a pole in the RHP for / > 0, which for 
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III. VELOCITY-FORCE CHARACTERISTICS 
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/ \ is situated at Sp = pl'^^ For d > 2 the pole will only exist in the RHP for / > 
= ld-2(d - 2)/2. For / \ /e, the position of the pole is given by Sp = A{d){f - fcf/^'^'^^ 
for 2 < d < A (where A(d) is a complicated function of d), whereas for d > A one has 
Sp = [{d — A)/{d — 2)]{f — fc). This abrupt change of behaviour for Sp ai d = A will have 
interesting consequences. 

Evaluating the contribution to 0(t) from the pole leads us to the following results for the 
velocity-force characteristics: i) for d < 2 one has for / \ the relation v{f) ~ j2/(2-d). 
ii) for d > 2, the velocity is zero unless / exceeds fc - one then has for / \ /c the relation 
v{f) ~ (Sf)^^^^ where 6f = f - and 

2<d<A 



a{d) = { (9) 
1 d>A, 

which signals the existence of an upper critical dimension du = A. 

In the absence of a pole in the RHP, the asymptotics of (j){t) arise from the contribution 
from the branch cut. The results for this case (along with the v{f) characteristics at the 
critical dimensions 2 and 4) are displayed for completeness in Table I. 

IV. SPATIAL PROFILES 

In this section we wish to evaluate the spatial profile of the interface for / > f^. We 
define the positive quantity Aif(x, t) = (A/2z/)(/i(0, t) — /i(x, t)). Making use of the Cole- 
Hopf transformation, along with Eq. (^, we have (after rescaling x ^ x = x/(4z/)) 

AiJ = - In ^ -i- + -f- ( dt' g{x, t' + Z^) 0(t - t') \ . (10) 



For f > fc and ft ^ 1, we have from the previous analysis (f){t) ~ e^*. In this case, the 
above expression reduces to 

Ai7 = -ln[/F(x)] (11) 

where 



LXJ 

F(x) = J dt' e-"*' (7(x, t' + f). (12) 



This expression is valid for |x| ^ v^/'^t. For this entire range we see that the height deviation 
A/i is independent of time. For extremely large distances from the forcing site (|x| ^ v^f^t), 
the influence of the force is insufficient to induce a velocity, and the surface distortion is 
essentially zero. 

We now concentrate on examining the form of the steady-state profile given by Eqs.(|TTl) 



and ([121) . In fact there are two spatial regimes for this profile, separated by the length scale 
Rc = f which will emerge as a critical radius. We shall denote the regime Rc ^ |x| ^ 
v^^'^t as region A, and the regime / ^ |x| ^ R^ as region B. It is important to stress that 
for vt ^ 1 and v <^ 1 (which is the case when / \ /c), regions A and B are both large and 
well-separated. Region A is less interesting and simpler to analyse, so we shall focus on this 
first. 

A. AH in region A 

In this region, -F(x) may be evaluated using a steepest-descents treatment. A straight- 
forward calculation yields the result 

I 2|x|/i?, + 0(ln(|x|/i?,)), d<2 
Aff(x) = i ' ' V VI 1/ c;;, ^^^^ 

[ 2|x|/i?, + O (ln(|x|/i?,)) + O (ln(i?,//)) , d>2 

The form of the corrections is of interest. We see that for d < 2, the leading term dominates 
for |x| ^ i?c as expected. This indicates that the overall height deviation within the critical 
radius is of order unity. For d > 2 the microscopic length scale enters into the corrections 
and the condition for the dominance of the leading term is now |x| ^ Rcln{Rc/l). We are 
therefore given the hint that for d > 2 the overall deviation of the height within the critical 
radius may be large (of order \n{Rc/l)). To see this exphcitly we now turn our attention to 
the form of AH for Ixl ^ Rc. 
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B. AH in region B 



Although the integral determining -F(x) appears quite simple, the analysis in region B 
is rather complicated. We therefore refer the reader to Appendix A where this analysis is 
described in some detail. The end result for AH in region B takes the form 

c{d)f |x|2-'^ + 0((|x|/i?,)2), d<2 
Ai/= <! 2/ln(|x|//) + 0(/), d = 2 (14) 

(rf-2)ln(|x|//) + 0(l), d>2 

where c{d) = 2T{d/2)/{2 — d). [It is instructive for d = 2 to rewrite the prefactor of the 
height deviation exclusively in terms of Rc - making use of the relation 2/ = 1/ ln{Rc/l)]. 

There are two main features of these results we wish to stress. First, there exists for all 
dimensions a critical radius Rc within which the surface profile undergoes a distortion from 
the linear form shown in Eq.(|r^). The critical radius scales with 5/ = / — /c as -Rc ~ (^/)" 
where u^d) = a{d)/2. In particular we note that for 2 < d < 4, u^d) = l/{d — 2); and for 
d> 4, u{d) = 1/2. 

The second important point is that by calculating the overall height distortion within 
the critical radius, i.e. AH{Rc), we see that 

0(1) 0< d<2 

AH{Rc) ~ " " (15) 

[ ln(i?c/0 d>2 
These results confirm the ideas propounded earlier, i.e. for d < 2 the distortion of the surface 
within the critical radius is negligible, whereas for d > 2 the distortion is a macroscopic 
rupture which scales logarithmically with distance. 

V. APPLICATIONS TO KPZ WEAK/STRONG-COUPLING TRANSITION 

As we mentioned in the Introduction, there are two main characteristics of the 
weak/strong coupling transition for the noisy KPZ equation (which arise from a RG study 
in powers of e = (i — 2) for > 2. These are the vanishing of the roughness exponent, which 
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implies a logarithmic roughness; and the existence of a correlation length ^ ~ ((? — gc)~^, 
where g is the dimensionless coupling constant, and u = l/e + 0{e)'^. There is a remarkable 
correspondence between these results, and those found in the previous section for the simple 
model of a deterministic KPZ interface under a local applied force. To reiterate these results, 
we found that for d > 2, the applied force induces a logarithmic distortion of the surface 
within a critical radius which scales as Rc ^ {f — fc)" with u = 1/e for 2 < d < 4, and 
u = 1/2 for d>4:. 

The natural question arises: are the critical properties of the noisy KPZ equation directly 
related to the simple problem discussed above? If so, then we have a much easier way 
to understand the weak/strong coupling phase transition, and also we may identify an 
upper critical dimension of 4 (for the unstable fixed point characterising the transition). A 
quantitative connection between the two sets of results is beyond our reach. However we 
shall make several points which certainly make the connection more plausible. Since we are 
interested in the weak/strong coupling transition, we henceforth restrict our attention to 
d>2. 

The coupling constant g in the noisy KPZ equation is proportional to the noise amplitude 
D which we shall take to be our control parameter. For D well below the critical value, the 
non-linearity in the KPZ equation is irrelevant. Thus, the surface behaves as an Edwards- 



Wilkinson (EW) interface |jTO|, which is smooth. As we increase the value of D, it becomes 
more and more likely for small regions in the interface to experience strong fluctuations 
due to rare events in the noise spectrum. If the smooth/rough transition is triggered by 
an inherent instability of the interface with respect to strong fluctuations, it is plausible 
to assume that such an instability is localized due to the local nature of rare events in the 
noise. Having established this, we may connect the two models by arguing that these strong 
local fluctuations resemble an effective force over long time scales. There are three essential 
steps in this argument which we discuss in turn. 

First, we must argue that the fluctuations appear as a constant force over some temporal 
interval r. Clearly any fluctuations in the 'downward' direction (i.e. anti-parallel to the 
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direction of positive h) will have no effect since they are completely suppressed by the non- 
linearity. We therefore consider a noise path (or noise 'history') which tends on average to be 
positive over the interval r within some region A of size l'^. This path may be considered as 
a force of magnitude / so long as it lies between the bounds / ± a. The fluctuations within 
these bounds are irrelevant. This fact may be established by considering the delta-function 
model Eq.@ but allowing for the force to be a stochastic variable. By power counting one 
may show that the fluctuations of the force are irrelevant (about the constant force flxed 
point) for d > 1. For the noise path to exceed / + a is extremely unlikely as the path is 
already a rare event even to have entered between these bounds. The path is much more 
likely to drop below f — c This is deemed as switching off the 'force'. The force will be 
switched on again when the noise path returns within the bounds. Thus the interval r over 
which the force acts is a random variable. 

Second, we must argue that this random switching off and on of the effective force is 
unimportant. We assume that the force has been able to create some disturbance in the 
interface over the region A, before it is switched off. How will the surface evolve? To answer 
this question one may simply solve the deterministic KPZ equation in the absence of any 
source, but with an initial condition corresponding to a localized disturbance. The result is 
that for any disturbance exceeding a critical height he ~ i^/A, the centre of the disturbance 
decays logarithmically slowly - it is essentially frozen. This effect is due to the strong 
influence of the non-linearity driving the disturbance upwards against the smoothing effect 
of the diffusion term. So during the periods of zero force, the disturbance in the region A is 
frozen and on switching on the force, the evolution continues as if there were no interruption. 

The third and flnal point is to argue that the effective force remains localized in the 
region A. This is in general not the case. The noise has short-range correlations in space 
and time and rare events are equally likely to occur anywhere in the surface. A simple 
resolution would be that once the rupture has been formed in A, then all other rare events 
in the vicinity may be rotationally averaged to give an effective force within A. A more subtle 
point is the following. Once the logarithmic distortion centered in A has been established, it 
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is essentially frozen during the periods of zero force. Therefore any rare event occurring in 
A will continue the evolution of the distortion. On the other hand, if the critical height for 
freezing is not small, then most rare events will fail to seed such a distortion of the interface, 
since the disturbances due to rare events which are of a height less than he will simply 
diffuse away. We therefore have a picture of extremely strong fluctuations being required to 
seed a distortion of height > he, with less strong fluctuations being required to sustain the 
evolution of the logarithmic rupture once it has been seeded. 

If this rough physical picture is correct then the correspondence between the local force 
model considered here, and the noisy KPZ equation at the weak/strong-coupling transition is 
estabhshed. The surface morphology resulting from such a picture is that of a dilute system 
of macroscopic distortions, each with a logarithmic profile. This is in contradistinction to the 
more conventional view of weak uniform fluctuations giving rise to logarithmic roughness. 
Such a difference in morphology should be discernible from numerical simulations of models 
believed to he in the KPZ universality class. 

[A physically plausible connection between these two models may be drawn also in the 
directed polymer representation. The strong-coupling phase of the KPZ model corresponds 
to the low-temperature phase of the directed polymer, in which the polymer is stongly local- 
ized onto low-energy paths, which have characteristic transverse fluctuations. Conversely, 
the weak-coupling phase of KPZ corresponds to the high-temperature phase of the directed 
polymer in which the polymer is liberated from the low energy path and performs thermal 
wandering. Hence the transition between these two phases for the directed polymer may 
be viewed as the localization transition to the lowest energy path (at this temperature). In 
this case the polymer is essentially subjected to a 'columnar' defect, although the defect 
is not of constant energy, and also has transverse fluctuations. Using similar arguments 
to those presented above in the interface language, one can motivate the idea that i) the 
energy fluctuations on the column are irrelevant within some bounds, and ii) the transverse 
fluctuations of the column are thermally averaged such that the effective column is straight 
on some transverse scale of order l'^. We believe a more quantitative connection may be 
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tractable in this directed polymer picture.] 



VI. CONCLUSIONS 

We have presented a non-perturbative analysis of the deterministic KPZ equation under 
the influence of a local driving force f, which for d > 2 must exceed a critical value fc in 
order to induce a non-zero velocity in the surface. The central result is that for d > 2 
and for / \ f^, the force creates a macroscopic distortion of the surface. This distortion 
has a logarithmic proflle and exists within a critical radius of size i?c ~ (/ ~ /c)" where 
z/((i) = 1/ {d—2) for 2 < d < 4 and = 1/2 for > 4. In the penultimate section we offered 
some physical arguments as to why one might consider this simple model to underlie the 
weak/strong-coupling transition in the noisy KPZ model. The main characteristics of such 
a scenario are that i) the logarithmic roughness at the transition is due to a dilute system of 
logarithmic distortions, as opposed to a uniform sea of fluctuations with logarithmic variance, 
and ii) the critical radius Rc is a physical realization of the correlation length found from 
the RG treatment. A more quantitative connection between these two models is presently 
beyond our reach, although one future possibility is to estimate the d-dependence of the 
weak/strong-coupling phase boundary based on the exact results obtained in this paper for 
the local force model, in conjunction with the statistical arguments presented above. 

It would be useful to apply the non-perturbative methods presented here to the case of 
a random local force. Such an analysis would allow one to tighten some of the arguments 
given in section 5. This problem is also of substantial interest in its own right as it may be 



mapped to a system of two directed polymers with random contact interactions |TT[| [12 
which has a wide range of applications. 
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APPENDIX A: 

In this appendix we analyse the function -F(x) in region B, i.e. for I <^ \x.\ <^ R^, = f "^Z^, 
in order to derive the relations shown in Eg. (p^ . The function F{x.) has the explicit form 



(cf. Eq.(0)) 

oo 

F(x) = j dt{t + f)-^/^ exp(-t;t) exp[-xV(t + 1% (Al) 



which may be re-expressed as 

F(x) = x^-'^ e^'' [h{x) + h{x)] (A2) 

where 

1 

h{x)= j dt r''/^ e-'^' e-^/\ (A3) 

and 

oo 

I^[X) = fdt g-at g-lA (^4) 



1 



where we have defined a = vx"^ = (x/Rc)"^ <^ 1 for notational convenience. 
The first integral may be simply evaluated to leading order in a. We have 

oo 

h{x) = fdu m'^/^-s g-« + 0(a) + 0(e-("/')'). (A5) 



The leading terms of the second integral may be extracted as follows. Referring to 
Eg. ([A^) , we expand the last factor of the integrand as a power series to give 

oo / -| \n 

U^) = E -^Jn{<^) (A6) 

where 

oo 

J„(a) = a^-^+^^Z" j dt (A7) 

a 

Integrating by parts yields the recursion relation 
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For d < 2, Jo(a) has the expansion 

Jo(a) = r(l - d/2) - (1 - d/2)-^ + 0(a) (A9) 

where T{z) is the Gamma function |]^. This form of Jq imphes Ji(a) = 2/d + 0(a'^/^) and 
for n > 1, J„ = (ra — 1 + d/2)^^ + 0(a). Returning to Eq. (|A6|) we have 



h{x) = r(l - t;/2) a'^/'-^ + C + 0(a'^/2) (AlO) 

where 

oo (_'\\n \ 
n=0 

Combining Eqs. (|A5|) and (|A10|) then gives for d < 2 

Ji + J2 = r(l - d/2) a^/2-1 _ (1 _ d/2)-^V{d/2) + 0(a'=^/2)_ (A12) 



From Table I we have for < 2 , = [/ 1(1 - d/2)Yl^'^-'^\ which in addition to Eqs.(^ 
and ( |A12| ) gives 



/F(x) = 1 - /(I - d/2)-^Y{d/2) |x|2"'^ + O ((|x|/i?e)') . (A13) 

Combining the above equation with Eq.(|Tl]) then yields the first relation in Eq.(|l^). 

A similar analysis suffices for higher dimensions, which we shall briefly sketch. For d = 2 
one finds 

+ /2 = -ln(a) -27 + 0(aln(a)), (AM) 

where 7 = 0.57721... is Euler's constant. From Table I we have for d = 2 , f = e~^^l^ ^ 
which in addition to Eqs.( [A2D and ( |A14| ) gives 

/F(x) = l-2/ln(|x|/0 + O(/). (A15) 

Combining the above equation with Eq.(|lTD then yields the second relation in Eq. ([I^) . 
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Finally for d > 2 one finds 



h + l2 = T{d/2 - 1) + O(a^) 



(A16) 



where f5 = d/2 ~ 1 for 2 < d < 4, and (3 = 1 for d > 4. From Table I we have for d > 2, 
f ^ f^ = [d/2 - ly-^, which in addition to Eqs.(^) and (|TOD gives 



/F(x) = r(d/2) (//|x|)'^-2 + O ((|x|/i?,)2'3(//|x|)^-2) . 



(A17) 



Combining the above equation with Eq. fpUD then yields the third relation in Eq . (|T4|) . 
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Table I 



d 


fc 




v{f) 


/</c 


f = fc 


f>fc 


0<d<2 





In (/V(2-<i)t) 





vt 


ci(d) /^/(^-'^) 


2 





-In (ln(t//2)) 





vt 


^exp(-l//) 


2 < d < 4 


(d-2)U-2 
2 ' 




^in(t//2) 


vt 




4 






In r 1 


vt 


1 (-J/Z/e) 

?^ ln(/c/<5/) 


d > 4 


2 *■ 




ln(t//2) 


vt 


(5///c) 



Caption for Table I 



Asymptotic behaviour of H{t) - the interface height at the forcing site - as a function of 
dimension d and applied force /. The symbol Sf is the excess force: df = f — fc- The 
constants Ci and C2 have the form ci(d) = [r(l - d/2)]2/(2-<i) and C2(d) = [r(2 - d/2)]2/(2-<i). 
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